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^ , Abstract 

pq , The result of the 3-year neutrino magnetic moment measurement at the Kalinin Nu- 

clear Power Plant (KNPP) with the GEMMA spectrometer is presented. Antineutri- 
no-electron scattering is investigated. A high-purity germanium detector of 1.5 kg 

f-» placed at a distance of 13.9 m from the centre of the 3 GW t h reactor core is used 

in the spectrometer. The antineutrino flux is 2.7xl0 13 P e /cm 2 /s. The differential 
method is used to extract v-e electromagnetic scattering events. The scattered elec- 
tron spectra taken in 5184+6798 and 1853+1021 hours during the reactor ON and 
OFF periods respectively are compared. The upper limits for the neutrino mag- 
netic moment [i v with and without atomic ionization mechanism were found to be 
5.0xl(r 12 /i B and 3.2xKT n /iB at 90% CL, respectively. 
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1 Introduction 



The Minimally Extended Standard Model predicts a very small magnetic mo- 
ment for the massive neutrino (//„ ~ lO" 20 ^) which cannot be observed in 
any experiment at present. On the other hand, there is a number of exten- 
sions of the theory beyond the Minimal Standard Model where the Majorana 
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neutrino magnetic moment (NMM) could be at a level of 10~( 10 - 12 ) /xb irre- 
spectively of the neutrino mass [1,2,3,4,5,6,7,8]. At the same time, from general 
considerations [9,10] it follows that the Dirac NMM cannot exceed 10 _14 /xb- 
Therefore, observation of the NMM value higher than 10~ 14 /xb would be ev- 
idence for New Physics and, in addition, would indicate [11,6,7,8] doubtless 
that the neutrino is a Majorana particle. 

It is rather important to make laboratory NMM measurements sensitive enough 
to reach the ~10 _11 /xb region. The Savanna River experiment by Reines' group 
could be considered the beginning of such measurements. Over a period of 
thirty years the sensitivity of reactor experiments only increased by a factor 
of three: from (2 . . . 4) x 1(T 10 /ib [12,13] to (6 . . . 7) x Kr n /XB [14,15]. Similar 
limits were obtained for solar neutrinos [16, 17], but due to oscillations at long 
distance (as well as matter-enhanced oscillations in the Sun) their flavor com- 
position changes and therefore the solar NMM results could differ from the 
reactor ones. 

In this paper, the results of a 3-year NMM measurement being performed 
by the collaboration of ITEP (Moscow) and JINR (Dubna) are presented. 
Measurements are carried out with the GEMMA spectrometer[18,19,15] at 
the 3 GWth reactor of the Kalinin Nuclear Power Plant (KNPP). 



2 Experimental approach 



A laboratory measurement of the NMM is based on its contribution to the v-e 
scattering. For nonzero NMM the v-e differential cross section is given [13] by 
a sum of the weak interaction cross section (dow/dT) and the electromagnetic 
cross section (do 'em / 'dT)\ 
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Fig. 1. Weak (W) and electromagnetic (EM) cross sections calculated for several 
NMM values. 
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where E v is the incident neutrino energy, T is the electron recoil energy, 8w 
is the Weinberg angle and tq is the electron radius (jrrl 



2.495 x 10" 25 cm 2 ! 



Figure 1 shows differential cross sections (1) and (2) averaged over the typical 
antineutrino reactor spectrum vs the electron recoil energy. One can see that at 
a low recoil energy (T -C E v ) the value of daw/dT becomes almost constant, 
while daEM/dT behaves as T _1 , so that the lowering of the detector threshold 
leads to a considerable increase in the NMM effect with respect to the weak 
unremovable contribution. 



To realize this useful feature in our GEMMA spectrometer [15], we use a 1.5 
kg HPGe detector with an energy threshold as low as 3 keV. To be sure that 
there is no efficiency cut at this energy, the "hard" trigger threshold is made 
twice lower (1.5 keV). 
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Fig. 2. Ge detector inside the active (Nal, PS) and passive (Cu, Pb) shielding. 

The background is suppressed in several steps. First, the detector is placed in- 
side a cup-like Nal crystal with 14 cm thick walls surrounded with 5 cm of elec- 
trolytic copper and 15 cm of lead (Fig. 2). Active and passive shielding reduces 
the external 7-background in the ROI to a level of ~ 2 counts/keV/kg/day. 
Being located just under reactor #2 of the KNPP (at a distance of 13.9 m 



from the reactor core, which corresponds to an antineutrino flux of 2.7 x 
10 13 z/ c /cm 2 /s), the detector is well shielded against the hadronic component of 
cosmic rays by the reactor body and technological equipment (overburden~70 
m w.e.). The muon component is also reduced by a factor of ~10 at ±20° with 
respect to the vertical direction and of ~3 at 70°. . . 80°, but a part of the resid- 
ual muons are captured in the massive shielding and thus produce neutrons 
which scatter elastically in Ge and give rise to a low-energy background. To 
suppress this background, the spectrometer is covered with additional plastic 
scintillator plates (PS) which produce relatively long /i-veto signals. Special 
care is taken to reduce non-physical low-amplitude circuit noise (afterpulses, 
radio frequency interference, microphonism, etc.). Thus, for example, the de- 
tector signal is processed by three parallel independent electronic channels 
with different shaping time (2, 4 and 12 /xs), which allows a primitive Fourier 
analysis [20] to be performed a posteriori, so artefact signals are discriminated 
(see the next section). 



3 Data taking and processing 



In order to get a recoil electron spectrum, we use a differential method com- 
paring the spectra measured during the reactor operation (ON) and shutdown 
(OFF) periods. In our previous work we considered Phase-I (13 months' mea- 
surement from 08.2005 to 09.2006, including 5184 and 1853 hours of the reactor 
ON and OFF periods, respectively). Today we can add Phase-II - 19 months 
from 09.2006 to 05.2008. Unfortunately, for some organizational and techni- 
cal reasons, there were several interruptions in the measurement. After the 
preliminary selection described below, 6798 ON-hours and 1021 OFF-hours of 
active time were found to be available for analysis. 

During the measurements, the signals of the HPGe detector, anticompton Nal 
shielding and outer anticosmic plastic counters, as well as the dead-time infor- 
mation, are collected on an event-by-event basis. The neutrino flux monitoring 
in the ON period is carried uot via the reactor thermal power measured with 
an accuracy of 0.7%. 

The collected data are processed in several steps. The first step involves dif- 
ferent selections aimed at suppressing nonphysical and physical backgrounds: 

(1) Bad run rejection. We reject those one-hour runs which correspond to 
the periods of liquid nitrogen filling and any mechanical or electrical work 
at the detector site, as they could produce noise. 

(2) Radioactive noble gas rejection. Unfortunately, the detector shield- 
ing turned out to be not tight enough against radioactive noble gases. 
To smooth away this design defect, we analyze energy spectra measured 



during each several hours and check the stability of the 7-background. If 
any visible excess of 81 keV ( 133 Xe), 250 keV ( 135 Xej or 1294 keV ( 41 Ar) 
7-line occurs, the corresponding runs are removed] 2 1 

(3) Detector noise rejection. For some obscure reasons our Ge detector 
happened to become noisy from time to time. In order to reject these 
noisy periods, the low-amplitude count rate is checked second by second, 
and those seconds which contain more than 5 events with E >2 keV are 
rejected. 

(4) Audio-frequency rejection. We reject those events which are separated 
by a time interval shorter than 80 ms or equal to (n • 20.0 ± 0.1) ms. 
In such a way we suppress the noise caused by mechanical vibrations 
( "ringing" ) and the 50 Hz power-line frequency. 

(5) Fourier rejection. As it was mentioned, the real and artefact signals have 
different Fourier spectra. To exploit this difference, we build three plots 
similar to that shown in Fig. 3: (E 2 vs Ei), (E 3 vs E 2 ) and (Ei vs £"3). 
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Fig. 3. Example of the Fourier analysis made with different shaping-times: ADC-2 
operates with 4 fis pulses, and ADC-3 operates with 12 fj,s pulses. Plot (a) is made 
before and (b) - after the "audio-frequency" rejection; it is seen that most of the 
rejected events are non-diagonal. (The color intensity scale is logarithmic.) 

The real signal falls into diagonals (Ei ~ E 2 — E 3 ) within the energy 
resolution, whereas any nonphysical artefact shows a different pattern. 
We only select diagonal events and thus additionally reject low- and high- 
frequency noise. To ensure the lowest cut-off, we replace Ei, E 2 and E 3 
by their linear combination E: 

E = aE 1 + bE 2 + cE 3 , (3) 

where the amplitudes a, b, c are chosen (subjectively) so as to make the 



2 In fact, these files are used later for the "noble gas" correction of the rest of the 
data. 
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vector E be antiparallel to the noise gradient (Fig. 3b). 

After the above rejections we construct energy spectra for the ON and OFF 
periods and correciLH them in two steps: 

(1) Noble gas correction. As our spectrometer is not located in a special 
laboratory, but in the technological room, there are sometimes short op- 
erational periods when the concentration of 41 Ar, 133 Xe and 135 Xe in this 
room becomes higher than usual. Spectra measured under these condi- 
tions are used to evaluate the contribution of each radioactive gas to 
the low-energy part of the background. These contributions normalized 
to the intensities of the corresponding 7-lines are then subtracted from 
those few ON and OFF spectra where small traces of these lines are still 
present. In this case the value of such correction in the ROo does not 
exceed 1-2%. 

(2) Low-energy threshold correction. The detection efficiency 77 just above 
the threshold E is measured with a pulser and fitted with the function 
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where a stands for the detector energy resolution. Experimental spec- 
tra are then corrected by the function (4) which in our case becomes 
significant at energies below 2.8 keV. 

As a result, we obtain energy spectra S for the ON and OFF periods which 
must be normalized to the corresponding active times Tqn and T OFF and then 
compared to each other, taking into account the additional neutrino-dependent 
term: 



^ = ^ + m^AW + X-EM) 



ON J- OFF 



(5) 



(effect) (background) (neutrino contribution) 



This term includes the fiducial detector mass rrid and the antineutrino flux <!>„ 
(known with an accuracy of 1.7% and 3.5%, respectively) multiplied by the 
sum of two neutrino contributions: the weak one (W), which can be calcu- 
lated easily by formula (1) and is completely negligible in our case, and the 



3 The corrections do not introduce a significant error in the final result as they 
affect both ON and OFF spectra in the same way. 

4 The Region-Of-Interest (ROI) in our analysis includes two fragments from 3 to 
9.4 and from 11.2 to 55 keV, i.e. the low-energy part of the continuous spectrum 
without peaks which could depend on the reactor operation. 



electromagnetic one (EM), which is proportional to the squared NMM value: 

Hi5%) 2 - (6) 



Unfortunately, the ON and OFF periods are not equal from the point of view 
of statistics. A usual OFF period is much shorter, and, therefore, the final 
sensitivity is limited by the background uncertainties. On the other hand, 
today, after three years of data taking, we know the ROI background structure 
with more confidence. It gives us the right to introduce additional information 
in our analysis, namely, to state that our background is a smooth curve. 

To implement this conventional idea, we fit the background OFF spectrum 
in the ROI from 2.9 keV to 55 keV with a parametrized smooth function (a 
sum of Gaussian, exponential and linear functions) or with splines; these fits 
produce slightly different results, and their spread is taken into account by 
the final systematic error. 

Then we compare the ON spectrum channel by channel with the obtained 
background curve and extract the X- value (or its upper limit) from Eq.(5). 
This evaluation is more complicated than expected because it is very difficult 
to count active times Ton and Toff precisely in a proper way (especially, after 
numerous selections of the events). To avoid possible errors caused by this 
procedure, we divide the active time normalization into two parts: absolute 
(Tqn) and relative (r = T n/T ff)- 

Roughly both the Tqn and Toff active times are estimated using several 
background 7-lines: the 238 keV line of 212 Pb, the 1173 keV and 1333 keV 
lines of 60 Co and the 1461 keV line of 40 K. This radiation originates from the 
pollution of the internal parts of the spectrometer and is therefore more-or-less 
stable in time and does not depend on the reactor operation. Comparing the 
intensities of the above lines measured with and without any selections, we get 
the estimates Tqn, Tqff with accuracy of 0.9% and 1.9%, respectively. This is 
not enough, however, to evaluate the r value with the required precision. We 
resolve this problem in the following way. 

The relative ON/OFF time factor r is represented as a product of its estimate 
f = Ton /Toff (which is a constant known with a bad accuracy of 2.1%) and a 
correction factor K (which should be not far from 1.0): r = Kr. Then Eq. (5) 
can be transformed to 

X ■ EM + W = (Son - K ■ f ■ S ff) • (T ON m d ^f 1 (7) 



As it is seen, the absolute time normalization Ton contributes to the final 
result in the same way as <& v or m d (i.e., simply as a factor) and therefore 



we can replace Ton by its estimate Ton- Standard systematic deviation 5X 
caused by this factor is not significant: 

5X 5(T ON ■ m d ■ $„) _ 
X Ton -md-^v 



(0.9%) 2 + (1.7%) 2 + (3.5%) 2 ~ 4.0% 
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Fig. 4. Example of two-dimensional charts x 2 {X,K) built for some portions of the 
GEMMA data. 

Considering X and K in Eq.(7) as two unknown free parameters and varying 
them around the point (X=0; K—l), we calculate residual sums x 2 and thus 
build a two-dimensional chart x 2 (X, i^) similar to the one shown in Fig. 4. 
The appearance of a closed graph in the chart corresponds to the minimum of 
X 2 ■ the non-zero eccentricity of the ellipse-like curve indicates some non-zero 
correlation between X and K, but this correlation is not too high, so that 
both X and K must be fitted only together, whereas can be preliminarily 
estimated separately. 

Indeed, according to Fig. 1, only the left-hand part of the ROI is really sensitive 
to the X-value, whereas the right-hand part is almost insensitive to X but 
still sensitive to the relative time factor r and hence to the K parameter. It 
means that preliminary time normalization can be performed not only with 
the background 7-lines, but also with a part of the continuous spectrum (e.g., 
from 20 to 55 keV). Both methods give very similar results, but the second 
one provides better precision. The differential ON-OFF spectrum built with 
this t estimate is shown in Fig. 5. 

The gain, energy scale, resolution and noise level were not exactly the same 
for different long-term ON measurement periods of Phase-II. As a result, the 
optimal coefficients a, b and c of Eq. (3), as well as the low-energy edge of 
the ROI, were also changed. The data of these periods must be analyzed 
separately and only then their results be combined in the final distribution. 
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Fig. 5. Fragments of the experimental ON and ON-OFF spectra measured in 
Phase-II (only a part of the available statistics is presented) . 

The analysis of the j-period consists in channel-by-channel comparison of the 
ON and OFF spectra over the ROI (with allowance for the weak contribution) 
and building a % 2 (X, X) chart similar to the one shown in Fig. 4. Projection 
(not a section!) of the whole chart onto the X-axis gives us the Xj value 
with a standard statistical deviation dXj. Taking into account the standard 
systematic deviation SXj given by Eq.(8), we get Xj±dXj±5Xj and build the 
probability curve Pj(X) which already includes the X-uncertainty. Multiplying 
probability curves of all analyzed periods (and also of both Phases I and II), we 
get the probability distribution P(X) = Ylpj(X) for the total measurement. 




Fig. 6. Extraction of the NMM limit with the required CL from the renormalized 
part of the probability distribution P(X). 

In order to estimate the part of systematic error introduced by the subjective 
choice of the background approximation function, (a,&,c)-coefficients and ROI 
edges, we repeat the above procedure several times for different reasonable 
combinations of these parameters and build a sum of the obtained distribu- 
tions (Fig. 6) which turns out to be about 7% wider. After a conventional 
renormalization recommended by the Particle Data Group [21] and described 
in our previous work [19], we extract the upper limit for the X parameter and 
thus get the following NMM limit: 



fi v < 3.2 x 10 _1 Vb (90%CL) 



(9) 



4 Atomic ionization effect 



When this work was almost ready, H.N. Wong et al.[22] noticed one more 
channel of the NMM interaction, namely, atomic ionization (AI) 
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The AI channel (Fig. 7) dominates at very low transferred energies T when 
the wave length of the virtual photon A (7*) = hc/T is comparable with the 
atomic dimension. In this case (i.e., at T ~ 1 keV), the virtual photon may 
interact not only with an atomic electron considered free, but also with the 
atom as a whole. According to [22] , the cross section of this process could be 
described as 



fit 



El 



Tta \i^bJ J- 
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where <j 1 a stands for the photoelectric cross section of a physical photon with 
the energy E~ = T. 





Fig. 7. The NMM interaction via the atomic ionization (AI) channel. Left: scheme of 
the process. Right: expected GEMMA count rates calculated for [i v = 3.2 x 10 _11 /iB 
and compared to the weak (W) and electromagnetic (u, e) scattering on free electrons 

(FE). 

The right-hand part of Fig. 7 shows the GEMMA count rates caused by weak 
(W) and electromagnetic (EMai and EM FS ) interactions calculated for \i v = 
3.2 x lCr n /iB- As one can see, the EMai is two to three orders of magnitude 
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highej 5 1 than the EMfe in the low-energy region. It means that from the 
point of view of the NMM-sensitivity, the region (T <10 keV) is especially 
important, at least while the wave length A (7*) is rather short with respect 
to the crystal structure. 
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Fig. 8. Differential ON-OFF spectra (points) measured in Phases I, II and expected 
count rates (curves) calculated with formula (11) for different NMM values in terms 
of KHVb- 

Unfortunately, our GEMMA spectrometer was not adjusted for the detection 
of T < 3 keV. Nevertheless, we estimated its sensitivity to the AI process. 
Figure 8 represents our experimental points measured in Phases I and II com- 
pared to the theoretical curves calculated with formula (11) for different NMM 
values. 
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Fig. 9. Probability distributions obtained for different subjective parameters used 
in the analysis; Ps represents their weighted sum. 

Exact mathematical analysis of both Phase I and Phase II results shows that 
the statistical precision of our data allows reaching (2 ... 3) x lCT 12 /^ (at the 
90%CL). At the same time, the possible systematic error (and, therefore, the 
wide spread of the separate curves of Fig. 9) is much higher. It originates 
mainly from the fact that the " hard" energy threshold of the detection system 



5 It should be mentioned that according to another approach [23], the AI effect is 
not so high and even may be negligible. 
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in Phases I and II was not low and stable enough. Under this condition, it is 
only possible to set the following NMM limit: 

\i v < 5.0 x KT 12 ^ (90%CL) (12) 



5 Conclusion 



The experimental NMM search with the GEMMA spectrometer has been going 
on at the Kalinin Nuclear Power Plant (Russia) since 2005. The HPGe detector 
of 1.5 kg placed 13.9 m under the core of the 3 GW t h water-moderated reactor 
is exposed to the antineutrino flux of 2.7 x 10 13 z/ e /cm 2 /s. As a result of the 
3-year measurement (about 13000 ON-hours and 3000 OFF-hours of active 
time), the upper limits of 5.0x10 -12 /ib and 3.2xlO _11 yUB at 90% CL were 
found for the NMM with and without using the Al-interaction mechanism, 
respectively. As it was mentioned by M.B. Voloshin [23] contrary to H.T. Wong 
et al. [22], the atomic contribution is not essential at all, so that the first of 
the above limits should be considered as illustrative only. 

At present, analysis of the data taken under improved conditions (Phase III) 
has just started, but it indicates that further sensitivity improvement of the 
spectrometer can be reached only by its significant upgrading. Within the 
framework of this new project (GEMMA-2) we will use the antineutrino flux 
of ~ 5.4 x 10 13 z/ e /cm 2 /s, increase the mass of the germanium detector by a 
factor of four and decrease the level of the background. The main improvement 
is expected to be the significant lowering of the energy threshold (below 1 ke V) . 
These measures will provide the possibility of achieving the NMM limit at the 
level of 1.0 x lO" 12 ^- 
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